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Abstract 

In the Standard Model, the charged current of the weak interaction is governed by a unitary quark 
mixing matrix that also leads to CP violation. Measurement of the Cabibbo-Kobayashi-Maskawa 
(CKM) matrix elements is essential to searches for new physics, either through the structure of 
the CKM matrix, or a departure from unitarity. We determine the CKM matrix element \V c b\ 
using a sample of 3 x 10 6 BB events in the CLEO detector at the Cornell Electron Storage Ring. 
We determine the yield of reconstructed B° — > D* + iu and B — > D*°£i? decays as a function of 
w, the boost of the D* in the B rest frame, and from this we obtain the differential decay rate 
dT/dw. By extrapolating dT/dw to w = 1, the kinematic end point at which the D* is at rest 
relative to the B, we extract the product |Vrf,|.F(l), where J-(V) is the form factor at w = 1. We 
find = 0.0431 ± 0.0013(stat.) ± 0.0018(syst.). We combine 1^1^(1) with theoretical 

results for .F(l) to determine \V cb \ = 0.0469 ±0.0014(stat.) ±0.0020(syst.) ±0.0018(theo.). We also 
integrate the differential decay rate over w to obtain B(B° -> D* + £u) = (6.09 ± 0.19 ± 0.40)% and 
B(B- -► D*°£u) = (6.50 ± 0.20 ± 0.43)%. 
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I. INTRODUCTION 



The elements of the Cabibbo-Kobayashi-Maskawa (CKM) quark mixing matrix (T|, |2[ are 
fundamental parameters of the Standard Model and must be determined experimentally. 
Measurement of the matrix elements tests unified theories that predict the values of these 
elements. It also offers a means of searching for physics beyond the Standard Model by 
testing for apparent deviations of the matrix from unitarity, deviations that could arise if 
new physics affected the measurement of one of its elements. The status of this test is often 
displayed using the famous "Unitarity Triangle" [0] . The CKM matrix element \V c b\ sets the 
length of the base of this triangle, and it scales the constraint imposed by tx (This constraint 
scales as |Vtf,| 4 ), the parameter that quantifies CP violation in the mixing of neutral kaons 

Two strategies are available for precise measurement of \V^\, both of which rely on the 
underlying quark decay b — > civ, where I indicates e~ or fi~. The first method combines 
measurements of the inclusive semileptonic branching fraction and lifetime to determine 
the semileptonic decay rate of the B meson, which is proportional to |K&| 2 - Theoretical 
quark- level calculations give the proportionality constant, thereby determining |V^,|, with 
some uncertainties from hadronic effects. This first approach relies on the validity of quark- 
hadron duality, the assumption that this inclusive sum is insensitive to the details of the 
various final states that contribute. 

The second approach uses the specific decay mode B — > D*£u or B — > Dip. The rate 
for these decays depends not only on |V^,| and well-known weak decay physics, but also 
on strong interaction effects, which are parameterized by form factors. In general, these 
effects are notoriously difficult to quantify, but because the b and c quark are both massive 
compared to the scale of hadronic physics, A « 0.5 GeV, heavy-quark symmetry relations 
can be applied to B — > D^Hu decays || |6|, 0, § H . In the limit mb, m c — > oo, the form factor 
is unity at zero recoil, the kinematic point at which the final state D^*> is a rest with respect 
to the initial B meson. Corrections to the infinite-mass limit are then calculated using an 
expansion in powers of A/m,Q. Luke showed |7j that the first-order correction vanishes for 
pseudoscalar-to-vector transitions, making D*lv decays more attractive theoretically than 
Div for \V ch \ determination. 1 Heavy Quark Effective Theory (HQET) JTD|, [TTj, 0, 0, [TJ 
exploits the heavy-quark symmetry and offers a rigorous framework for quantifying the 



hadronic effects with relatively small uncertainty [15, 16]. 

In this paper, we report more fully on a recently published |T7] measurement of |V^{,| using 
B — > D*iv decays that are detected in the CLEO II detector at the Cornell Electron Storage 
Ring (CESR). The B — > D*iv decays are fully reconstructed, apart from the neutrino. The 
analysis takes advantage of the kinematic constraints available at the T(4S) resonance, 
where the data were collected, to suppress backgrounds, help distinguish D*iu from similar 
modes such as Diiu, and provide superb resolution on the decay kinematics. This analysis 
is the first since a previous CLEO result |18| to use not only B° — > D* + lv decays, but also 



B~ — > D*°iu decays |T{|. Consistency between these two modes is a valuable cross-check of 
our results. 

We reconstruct D* + candidates and their charge conjugates (charge conjugates are im- 



There are experimental advantages as well: a larger branching fraction, a distinctive final state with 
the narrow D*, and less phase-space suppression than the P-wave decay B — > Dlv near the important 
zero- recoil point. 
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plied throughout this paper) through the modes D* + — > D°tt + and D° — > K~tt + , and we 
reconstruct D*° candidates through the modes D*° — > D tt°, D° — > A" _ 7r + , and 7r° — > 77. 
Each Z)* candidate is combined with an electron or muon candidate. We then divide the 
reconstructed candidates into bins of w, where w is the scalar product of the B and D* 
four-velocities, and equals the relativistic 7 of the D* in the B rest frame. 2 Given these 
yields as a function of w, we fit simultaneously for parameters describing the form factor 
J~{w) and the normalization at w = 1. This normalization is proportional to the product 
|Vd>| 2 .F 2 (l), and combined with the theoretical results for J 7 (I), it gives us \V c b\. 



II. EVENT SAMPLES 

Our analysis uses 3.33 x 10 6 BB events (3.1 fb" 1 ) produced on the T(4S) resonance 
at the Cornell Electron Storage Ring (CESR) and detected in the CLEO II detector. In 
addition, the analysis uses a sample of 1.6 fb -1 of data collected slightly below the T(45) 
resonance for the purpose of subtracting continuum backgrounds. Because of miscalibration 
of low-energy showers in the calorimeter in a subset of the data, we use only 3.04 x 10 6 BB 
events (2.9 fb -1 ) produced on the T(4S') resonance and 1.5 fb -1 of data collected below the 
T(4S*) resonance for reconstructing B~ — > D*°£v candidates. 



The CLEO II detector [g(J has three central tracking chambers, immersed in a 1.5 T 
magnetic field, that measure charged particle trajectories and momenta. The momentum 
resolution is 5 MeV/c (12 MeV/c) for particles with a momentum of 1 GeV/c (2 GeV/c) 
(typical for the lepton and the K and 7r from the D°) and 3 MeV/c for particles with 
momentum less than 250 MeV/c (typical for the 7r + from the D* + ). A CsI(Tl) calorimeter 
surrounds both the tracking chambers and a time-of-flight system that is not used for this 
analysis. The calorimeter provides photon detection and assists with electron identification. 
The energy resolution of the calorimeter is 3.8 MeV for 100 MeV photons, a typical energy 
for photons from the decay of the 7r° from the D*° decay. The outermost detector component 
consists of plastic streamer counters layered between iron plates and provides detection of 
muons. 



We also use simulated event samples from a GEANT-based Monte Carlo simulation. 
With this Monte Carlo, we produce large samples of simulated B — > D*iv decays as well as 
a sample of 16 x 10 6 BB events to study some backgrounds. 



III. EVENT RECONSTRUCTION 



The T(4S') is produced and decays at rest, and each daughter B meson is produced with 
a momentum of about 0.3 GeV/c. As a result, T(45) — > BB events tend to be isotropic, 
or "spherical," with particles carrying energy in all directions. When the electron-positron 
collisions in CESR do not produce T(4S')'s, they can produce, among other things, qq quark 
pairs, where the q is a c, s, u, or d quark. Because the mass of these quark pairs is much 
lower than the energy of the beam, the daughter particles of these quarks' hadronization have 
higher momenta than the B's. These events tend to have a more "jetty" appearance; that 



2 The variable w is linearly related to q 2 , the squared invariant mass of the virtual W, via w — {m 2 B + 
mfjt — q 2 ) / (2niBmD*) , where tub and rac> are the B- and D*-meson masses. 
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FIG. 1: The m(Kir) distribution for (a) D* + £u and (b) D* £v candidates. All requirements are 
met except \m(K-K) — 1.865| < 0.020 GeV/c 2 . We accept candidates that fall between the arrows. 



is, the energy in the event tends to be distributed back-to-back. The ratio of Fox- Wolfram 
moments H2/H0 measures an event's jettiness, with values of the ratio approaching zero 
for spherical events, and approaching one for jetty events. To suppress non-BB events, we 
require that the ratio of Fox- Wolfram moments H2/H0 be less than 0.4, a condition satisfied 
by 98% of BB events containing a D*£u decay. 

To reconstruct D*£v candidates we first form D° — > K~tt + candidates from all possible 
pairs of oppositely charged tracks, alternately assigning one the kaon mass and the other the 
pion mass. We require a fiducial cut of | cos(9| < 0.9 for tracks, where 9 is the polar angle of 
the track's momentum vector with respect to the e + e~ beam axis. Tracks outside this fiducial 
region are excluded from consideration because they are poorly measured, having passed 
through the endplate of one of the inner tracking chambers and therefore either traversing a 
significant amount of material before entering the outer tracking chamber or never entering 
it at all. We reconstruct the invariant mass m(Kir) of the D candidate with a resolution 
of about 7 MeV/c 2 , accepting candidates that lie in the window \i7i(Ktt) — 1.865| < 0.020 
GeV/c 2 . The m(Kir) distributions for D* + £v and D*°£u candidates are shown in Fig. [I]. 

The pions produced in the decay D* — > Dtt have low momentum (< 250 MeV/c) because 
the combined mass of the D° and tt is within 8 MeV/c 2 of the mass of the D* . We label 
these pions "slow." For D* + candidates, we add a slow 7r + candidate to a D° candidate, 
requiring that the slow pion have the same charge as the pion from the D decay. This pion 
must also satisfy | cos#| < 0.9. The K and 7r are fit to a common vertex, and then the slow 
7r + and D° are fit to a second vertex using the beam spot constraint. For this vertexing we 



use error matrices from our Kalman fitter 



We then form Am = m{Kin:) — m(K7r). 
We look at Am rather than m(K7nr) because subtracting the D° candidate mass from the 
D* candidate mass cancels some errors in reconstructing the D°. A plot of Am for D* + 
candidates is shown in Fig. 0(a). The vertex constraints improve the resolution by about 
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FIG. 2: The Am distribution for (a) D* + £u and (b) D*®tv candidates. All requirements are met 
except \Am - 0.14544| < 0.002 GeV/c 2 for D* + £v and \Am - 0.1422| < 0.003 GeV/c 2 for D*°£i?. 
We accept candidates that fall between the arrows. 



20% to 0.7 MeV/c 2 . We require |Am - 0.14544| < 0.002 GeV/c 2 for D* + candidates. 

For D*° candidates, we add a slow tt° — > 77 candidate to the D° candidate. We construct 
771(77) f° r slow 7T° candidates from showers in the Csl calorimeter whose position is incon- 
sistent with extrapolation of any of the tracks reconstructed in the event. We require that 
the lateral pattern of energy deposition in the calorimeter be consistent with expectations 
for a photon. Particles with |cos#| > 0.71 travel through the endplate of the outermost 
tracking chamber before reaching the calorimeter, again traversing a significant amount of 
material. We therefore require that both photon candidates satisfy | cos^l < 0.71 so as to 
remain in the part of the calorimeter with the best energy and position resolution. Both 
photons must have energy greater than 30 MeV to limit background from soft showers. We 
also require the invariant mass 777(77) to give the known tt° mass within roughly three times 
the resolution of 5 MeV/c 2 : 0.120 GeV/c 2 < 777(77) < 0.150 GeV/c 2 . The Am resolution 
for D* 0, s is about 0.9 MeV/c 2 , so we require |Am - 0.1422| < 0.003 GeV/c 2 . The Am 
distribution for D*° candidates is shown in Fig. 0(b), and the 771(77) distribution is shown 
in Fig. 0. 

Finally, we require the momentum of the D* candidate to be less than ^^E B — miKixTx) 2 , 
(approximately 2.5 GeV/c), where Eb is the energy of the beam. This requirement sup- 
presses background from non-BB events. 

We next combine the D* candidate with a lepton candidate, accepting both electrons 
and muons. Electrons are identified using the ratio of their energy deposition in the Csl 
calorimeter to the reconstructed track momentum, the shape of the shower in the calorimeter, 
and their specific ionization in the tracking chamber. We require our candidates to lie in 
the momentum range 0.8 GeV/c < p e < 2.4 GeV/c, where the upper bound is the end point 
of D*lv decays. This momentum selection is approximately 93% efficient for B — > D*e~v 
decays. We require muon candidates to penetrate two layers of steel in the solenoid return 
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FIG. 3: The 772(77) distribution for D*°£u candidates. All requirements are met except 0.120 
GeV/c 2 < 771(77) — 0.150 GeV/c 2 . We accept candidates that fall between the arrows. 



yoke, or about 5 interaction lengths. Only muons with momenta above about 1.4 GeV/c 
satisfy this requirement; we therefore demand that muon candidates lie in the momentum 
range 1.4 GeV/c < p e < 2.4 GeV/c. This more restrictive muon momentum requirement 
has an efficiency of approximately 61%. We require both muon and electron candidates 
to be in the central region of the detector (|cos#| < 0.71), where efficiencies and hadron 
misidentification rates are well understood. The charge of the lepton must match the charge 
of the kaon, and in the case of D* + £u decays, be opposite that of the slow pion. 

The remaining reconstruction relies on the kinematics of the B — > D*£v decay. We 
first reconstruct cos8b-d*£, the angle between the -D*-lepton combination and the B meson, 
computed assuming that the only unreconstructed particle is a neutrino. This variable helps 
distinguish B — > D*£v decays from background and is necessary for the reconstruction of w. 
To form cos 9 b-d*i, we first note that the 4-momenta of the particles involved in B — > D*£v 
decay are related by 

pI = (Pb~Pd* ~Vtf- (1) 
Setting the neutrino mass to zero gives 

= m 2 B + m{D*£) 2 - 2 [E{B)E{D*£) - p(B) ■ p{D*£)] . (2) 

We solve for the only unknown quantity, the angle between the B meson and the vD*-lepton 
pair: 

2E(B)E(D"e) - m| - m(D'tf 

cos9fl - D " = wmvm ■ (3) 

In forming cos8b-d*£, we use the momenta of the D* and lepton candidates as well as the 
B mass and average B momentum, measured in our data. At CESR, a symmetric e + e~ 



collider operating on the T(4S') resonance, the B energy and therefore momentum is given 
by the energy of the colliding beams. Instead of relying on beam energy measurements 
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FIG. 4: The difference between the reconstructed w and generated w for simulated D* + lv (a) 
and D*°£u (b) decays in the intervals 1.1 < w < 1.15 (solid) and 1.4 < w < 1.45 (dashed). The 
normalization of all four histograms is arbitrary. 



based on storage ring parameters and subject to significant uncertainties, we determine the 
average B momentum directly using fully reconstructed B decays to hadrons. The energy 
spread of the beams and run-to-run energy variations lead to a distribution of B energies 
and momenta. By measuring the momentum distribution of fully reconstructed hadronic B 
decays in our data sample, we determine the energy spread intrinsic to CESR, which is then 
used to simulate BB pair production in our Monte Carlo. For cos Ob-dh we use the true D* 
mass rather than the reconstructed m(K7nr) to avoid a bias in the cos8b-d*£ distribution 
of the Am sideband, which we use to determine a background. 

We next estimate w for each candidate. Exact reconstruction of w, the boost of the D* 
in the rest frame of the B, requires knowledge of the B momentum vector. Although the 
magnitude of the B momentum is known, the B direction is unknown. However, it must 
lie on a cone with opening angle 6b-d*i around the D*£ direction. We calculate w for all 
B flight directions on this cone and average the smallest and largest values to estimate w, 
with typical resolution of 0.03. We divide our sample into ten equal bins from 1.0 to 1.5, 
where the upper bound is just below the kinematic limit of 1.504. For a few candidates, the 
reconstructed w falls outside our range; we assign these to the first or last bin as appropriate. 
Figure ^ shows the distributions of reconstructed w minus generated w in the third and ninth 
w bins for simulated B° — > D* + lv and B~ — > D*°£u decays. 

In the high w bins, we suppress background with minor loss of signal efficiency by restrict- 
ing the cosine of the angle between the momenta of the D* and of the lepton (cos^d*-^). 
The distribution of cos9o*-i versus w is shown in Fig. ^ for simulated B — > D*lv decays. 
Some backgrounds are uniformly distributed in this angle. The accepted angles are listed in 
Table § 
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TABLE I: The accepted regions of the cosine of the angle between the D* and the lepton in each 
w bin. 



w bin 


w limits 


min. 


Accepted cos9d*-£ 

max. 


1-5 


< 1.25 


-1.00 


1.00 


6 


1.25-1.30 


-1.00 


0.25 


7 


1.30-1.35 


-1.00 


0.00 


8 


1.35-1.40 


-1.00 


-0.25 


9 


1.40-1.45 


-1.00 


-0.50 


10 


> 1.45 


-1.00 


-0.75 



IV. EXTRACTING THE D*lv YIELDS 
A. Method 

At this stage, our sample of candidates contains not only B — > D*£u decays, but also 
B — > D**iv and B — > D*-nlv decays and various backgrounds. In the following, we refer to 
B — > D**£u and non-resonant 5 — > D*n£u decays collectively as S -> D*Xlv decays. In 
order to disentangle the Z}*£z/ from the D*X£u decays, we use a binned maximum likelihood 
fit [25| to the cos8b-d*i distribution. As shown in Fig. || B — > D*lv decays are concentrated 



in the physical region, —1 < cos 9b-d*£ < 1> while the missing mass of the D*X£u decays 
allows them to populate cos 9b-d*i < — 1- In this fit, the normalizations of the various 
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FIG. 6: The cos9b-d*£ distributions for simulated D*£v and D*Xtv decays. 



TABLE II: The regions of cos 9 b-d*i over which we perform a binned maximum likelihood fit. 



w bin 


w limits 




cos8b-d*£ fit region 






min. 


max. 


1-6 


< 1.30 


-8.0 


1.5 


7 


1.30-1.35 


-6.0 


1.5 


8 


1.35-1.40 


-4.0 


1.5 


9 


1.40-1.45 


-3.0 


1.5 


10 


> 1.45 


-2.0 


1.5 











background distributions are fixed and we allow the normalizations of the D*£u and the 
D*X£u components to float. For each w bin, we fit over a cos8b-d*£ region chosen to 
include 95% of the D*X£v events in that bin. These regions are listed in Table |TJ. 

The distributions of the D*£v and D*X£v decays come from Monte Carlo simulation. 
We simulate D*£v decays using the form factor of [EB] and include the effect of final-state 
radiation (B — > D*iv^) using Photos P7| . For D*X£v, we model D**£v modes according 
to ISGW2 |gSJ and non- resonant D*ix£v from Goity and Roberts pi]. O ur model for D*X£v 
is dominated by approximately equal parts of D\£v and D*7r£v. The other backgrounds, 
and how we obtain their cos8b-d*£ distributions and normalizations, are described in the 
next section. 



B. Backgrounds 

There are several sources of decays other than B — > D*£v and B — > D*X£v. We divide 
these backgrounds into five classes: continuum, combinatoric, uncorrelated, correlated, and 
fake lepton. As an indication of the relative importance of the various backgrounds, in 
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TABLE III: The contribution of each background Bi to the total number of candidates 7V tot in the 
range — 1 < cos 6 b-d*i < 1 for the D* + £v and D*°£v analyses. The relative size to D*£u signal 
Bi/S is also given for the same cos 6 b-d*i interval. 



Background 


D* + £v Contribution 


D £u Contribution 




Si/Nut (%) 


Bi/S (%) 


Bi/N tot (%) 


Bi/S (%) 


Continuum 


3.8 


4.7 


2.8 


5.1 


Combinatoric 


7.9 


10 


38 


70 


Uncorrelated 


4.4 


5.6 


4.7 


8.6 


Correlated 


0.4 


0.5 


0.1 


0.2 


Fake Lepton 


0.5 


0.6 


0.2 


0.4 



Table |f| we give both the fraction Bi/N tot of candidates from each background source 
relative to all candidates and the ratio Bi/S of each background source to D*£v signal. 
Because signal events populate the physical region —1 < cos 9b-d*£ < 1; we compute both 
Bi/Ntot and Bi/S using only candidates in this "signal region." We discuss each background 
and how we determine it below. 



1. Continuum background 

At the T(45) we detect not only resonance events (T(4S ( ) — ► BB), but also non-resonant 
events such as e + e~ — > qq. This background contributes about 4% of the candidates within 
the signal region for D* + £u decays, and about 3% for D*°£u decays. This is about 5% 
relative to the D*£v signal. In order to subtract background from this source, CESR runs 
one-third of the time slightly below the T(4S) resonance. For this continuum background, we 
use the cos9b-d*i distribution of candidates in the off-resonance data scaled by the ratio of 
luminosities and corrected for the small difference in the e + e~ — > qq cross sections at the two 
center-of-mass energies. In reconstructing cos9b~d*£i we scale the energy and momentum 
of the D* and lepton by the ratio of the center-of-mass energies and use the B momentum 
measured in on-resonance data to compute the B energy. This continuum background 
includes combinatoric and fake lepton backgrounds arising from continuum processes. 



2. Combinatoric background 

Combinatoric background candidates are those in which one or more of the particles in 
the D* candidate does not come from a true D* decay. This background contributes 8% 
of the candidates in the signal region for D* + £u; for D*°£u, which suffers from random 
shower combinations and does not benefit from the charge correlation of the slow pion, this 
background contributes 38% of the candidates in the signal region. Relative to the D*£v 
signal, the combinatoric background is 10% for D* + £u and 70% for D*°£v. 

The cos 9 b-d*z distribution of the combinatoric background is provided by D*-lepton 
combinations in the high Am sideband. We choose the Am sidebands of 0.155 GeV/c 2 
< Am < 0.165 GeV/c 2 for D* + £v and 0.147 GeV/c 2 < Am < 0.165 GeV/c 2 for D*°£v. For 
values of Am above these ranges, the slow pions tend to be faster, and therefore cos9 B -dh 
tends to be larger, while regions closer to the Am signal region include signal decays in 
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FIG. 7: The Am distribution of candidates in the third w bin (1.1 < w < 1.15) for (a) D* + lv 
candidates and (b) D*°lv candidates with the result of the fit superimposed. The data (solid circles 
or squares) are superimposed with the combinatoric background distribution (dashed curve) and 
the sum of the background and the D* signal (solid histogram). In (a) the arrows delimit the fit 
region. In (b), the shaded histogram shows combinations in which only one of the two photons 
forming the 7r° candidate was correct. Unless indicated otherwise, the error bars provided in all 
figures are statistical only. 



which the slow pion is poorly reconstructed. With this choice, only 3.5% and 0.4% of the 
D*iv decays fall in the sideband for D* + £v and D*°£u, respectively. 

The normalization of the Am sideband candidates is determined in each w bin from a fit 
to the Am distribution with the sum of properly reconstructed D*'s and the combinatoric 
background. The line-shape for the D* peak is taken from simulated D*lv decays. The 
D*°£is line-shape includes D*° candidates in which only one of the two photons constituting 
the 7T° was correct. Since these candidates preferentially populate the Am signal region, a 
few (3.9% of all D*°£v decays) remain after our combinatoric background subtraction and 
are included in our D*®ii> signal. For D* + £u we assume a background shape of the form 

n(Am — m n ) a exp { [ci(Am — m n ) + c 2 (Am — m^) 2 ] } , (4) 

where c% and c 2 are constants fixed using an inclusive D* + sample, and we vary n, a, and 
the normalization of the signal peak. For D*®lv we assume a background shape of the form 

n(Am — m n ) a exp [b( Am — m T )] (5) 

and vary n, a, b, and the normalization of the signal peak. The fits for D* + £u and D* Q £i> 
are shown for a representative w bin in Fig. |7|. The normalizations are shown in Fig. [|. 

As a test of this background estimate, we carry out the same procedure used in data on 
a sample of 16 million simulated BB events. Because combinatoric background originates 
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FIG. 8: The ratio of the number of combinatoric background candidates in the Am signal region to 
the number in the Am sideband, as determined from fits to Am for (a) D* + lv and (b) D*°£u. The 
cos6b-d*£ distribution of the combinatoric background is provided by the sideband candidates 
normalized by this ratio. 



from random combinations of tracks and showers, we expect our Monte Carlo, which is 
tuned to reproduce track and shower multiplicity and momentum distributions of B decays, 
to provide a reliable check of the background estimation procedure. We compare the true 
background in the Am signal region with the background estimate formed using the Am 
sideband region. There is a concern that kinematic differences between candidates in the Am 
signal and sideband regions could cause a difference in the cos 8 ] b-d*i shape of the estimated 
and true backgrounds. Figures |S] and |TD| show the true and estimated backgrounds for the 
Monte Carlo sample. We observe that the shapes do differ for D* + £u, consistent with the 
effect of the strong momentum dependence of the slow-pion efficiency (Section |V1 B| ). The 
agreement is better for D*°£u. We evaluate the systematic error from this sideband technique 
in Section |V1A2 . 

This method of background estimation overlooks a small component of the combinatoric 
background, a component that arises from D* decays in which the slow pion is properly 
found but the D° candidate is constructed from the products of a D° decay other than 
D° — ► K~n + . Although the D° is misreconstructed, this background will still peak in the 
Am signal region. Most of these candidates have m(KiT) below our signal region, but Monte 
Carlo simulation shows that D° — > p + n~ / p~n + decays could contribute a few candidates to 
the m(K7r) and Am signal regions. Although these decay modes have not yet been observed, 
a combined branching fraction of 1.3% is plausible given the measured branching fraction of 
all D° — > 7t + 7t~tt° decays; with this branching fraction, these modes would increase our D*£v 
yield by (0.3 ± 0.2)%. In a sample of 16 million simulated BB decays, several other modes 
also contribute, bringing the total contribution to (0.5 ±0.3)%. The contributing modes are 
listed in Table [TV]. As this contribution has little effect on our results, and as the branching 
fractions of the main contributing modes (D° — > p + 7i~ and D° — > p~n + ) are unmeasured, 
we account for it in the combinatoric background systematic error, but otherwise neglect it. 
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TABLE IV: Decay modes of the D, other than D — ► K~tt + , that are not fully subtracted by the 
Am sideband. The third column shows the expected contribution (relative to D° — > K~tt + ) to 
the D*£u yield from each mode in the m(Kir) and Am signal regions. 



Mode 


Branching Fraction(%) 


Contribution(%) 


D° -> K + K~ 


0.425 ±0.016 a 


0.05 ± 0.03 


D° -> vr+vr- 


0.152 ±0.009 a 


0.02 ± 0.02 


£>° -> if-|it/ 


3.47 ±0.17 a 


0.02 ± 0.02 


£)0 _> tt+vt-ttO 


1.6 ± 1.1 ab 


0.33 ± 0.24 




10.8 ± 1.0 a 


0.02 ± 0.02 


D° -> A*-vr+ 


1.7 ±0.2 a 


0.01 ±0.01 




0.37 ±0.06 ac 


0.07 ±0.05 


Total 




0.52 ±0.25 



"From Ref . 

b The simulation includes nonresonant D° — > Tr + ir~ir° and resonant D° — > /J7r submodes. 
c Assuming lepton universality, we use the Z) — > 7r~e + ^ branching fraction for D — > ir~ tv. 



3. Uncorrelated background 

Uncorrelated background arises when the D* and lepton come from the decays of differ- 
ent B mesons in the same event. This background accounts for approximately 5% of the 
candidates in the signal region for both D* + £u and D*°£v decays, contributing 6% and 9% 
relative to D* + £v and D* Q £v, respectively. We obtain the cos 9 B _ D *i distribution of this 
background by simulating each of the various sources of uncorrelated D*'s and leptons and 
normalizing each one based on rates measured from or constrained by the data. We classify 
the D* and the lepton according to their respective sources because different sources give 
different momentum spectra for the D* and lepton, and therefore different distributions in 

cos 9 B -D*£- 

There are three components of uncorrelated background that contribute to both the 
D* + £u and D*°£u modes. The first component consists of a lower- vertex D* (i.e., from 
b — > c transitions) combined with a secondary lepton (i.e., from £> — c — ^ s£v) (primary 
leptons from the other B have the wrong charge correlation); this is the largest component 
for the D* + £u mode. Secondly, uncorrelated background can also occur when the B° and B° 
mix or when a D* from the upper- vertex (i.e., from b — > c, as in b — > ccs) is combined with 
a primary lepton (i.e., from b — > c[u)£v). Finally, in the D*°£v case, the largest source of 
uncorrelated background consists of candidates in which the K and 7r from a lower- vertex D* 
have been exchanged and paired with a primary lepton from the other B. (This background 
does not occur for D* + £v because we constrain the charge of the slow-pion candidate to be 
opposite to that of the kaon.) 

We first determine the production rate of upper-vertex _D*'s from B decays using the 
measured branching fractions of modes such as B — > DWjW^"" ||3i||. \y e do this in two D* 



momentum bins, relying on our simulation of such decays for the D* momentum distribution. 
To determine the lower-vertex D* production rate, we measure the rate of inclusive D* 
production from B decays in the data in each momentum bin and subtract the upper-vertex 
contribution from each. The results are shown in Table 0. We determine the background 
contribution from D's reconstructed with exchanged A°s and 7r's by studying inclusive D* + 
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TABLE V: The rate per BB pair used to normalize the D* elements of the uncorrelated background. 
The errors indicate the variation used to assess the systematic uncertainty in the background. 

Rate p D * < 1.3 GeV/c p D > > 1.3 GeV/c 

lower-vertex, D* + 0.281±0.032 0.242±0.015 

lower-vertex, D*° 0.231±0.031 0.272±0.014 

upper-vertex, D* + 0.048±0.024 0.012±0.006 

upper-vertex, D*° 0.048±0.024 0.004±0.002 



TABLE VI: Modes that contribute to the correlated background, their assumed branching fractions 
(BF), and fraction of the total correlated background. The numbers given are for the D* £i> mode. 



The contributions to the D* + H,u mode are similar. 

Mode BF || (%) Fraction(%) 

B -»• D*°Xt~p L65 36\0 

B -> D*~D*° 3.08 18.7 

B -> D*°DWkW 5.6 18.7 

B -> D*°Z>* 3.4 15.1 

5 -> DjD*° 1.39 6.5 

5 -> L>* 7 X; 7 -> e+e" 3.6 

5 -> D*°tt-; ir- -» ^t-p 0.47 1.4 



decays with the charge correlation of the slow pion reversed. 

We normalize the primary lepton decay rate for leptons with momenta between 0.8 GeV/ c 
and 2.4 GeV/c to its measured value of (8.99±0.42)% where the error includes statistical 
and systematic errors; since this measurement was made at CLEO, we include only the 
systematic errors that are uncorrelated with our analysis. Likewise, we adjust the secondary 
lepton rate for leptons with momenta between 0.8 GeV/c and 2.4 GeV/c to its measured 
value of (1.53 ± 0.12)% [^]. Finally, we adjust Xd, the B° — B° mixing probability, to its 
measured value of 0.174 ± 0.009 pOl . 



4- Correlated background 

Correlated background candidates are those in which the D* and lepton are decay prod- 
ucts of the same B, but the decay was not B — > D*tv or B — ► D*Xiv. The most common 

~ - ~ (*)— 

sources are B — > D*t~u followed by leptonic r decay, and B — > D*D S followed by semilep- 

tonic decay of the D~. The uncorrelated background contributes 0.5% and 0.2% compared 

to the D* + £u and D*°£u signals, respectively. The background is small; we therefore rely on 

our Monte Carlo simulation to quantify it. The decay modes and branching fractions used 



are listed in Table VI 
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5. Fake lepton background 

Fake lepton background arises when a hadron is misidentified as a lepton and is then 
used in our reconstruction. Fake leptons make up 0.5% of candidates in the signal region 
for D* + £v and 0.2% for D*°£u; relative to signal, the background contributions are about 
0.5%. To assess this background we repeat the analysis, using hadrons in place of the lepton 
candidates. After subtracting continuum and combinatoric backgrounds, we normalize the 
cos8b-d*i distributions with the probability for a hadron to fake an electron or muon. We 
measure the momentum- dependent fake probability using kinematically identified samples 
of hadrons: pions are identified using Kg — > 7r + 7r~ decays, kaons using D* + — > D°tt + — > 
K~7T + 7r + , and protons from A — > pn~ . The fake probabilities are then weighted by species 
abundance in B decays and the momentum spectrum of hadronic tracks in events with an 
identified D* + to obtain an average fake rate of 0.035% for a hadronic track to fake an 
electron and 0.68% to fake a muon. 



6. D*lu and D*Xlv cos Qb-d*i distributions 

The cos 9 b-d*£ distributions of D*£v and D*X£v decays are obtained from simulated 
BB events in which one of the B's is required to decay to D*£v or D*X£v. Since the other 
B in the event also decays, the cos 8b-d*£ distributions can contain the same backgrounds 
listed above. Using generator-level information, we veto all background sources except the 
combinatoric background, for which we perform the same Am sideband subtraction used in 
the data. In the sideband subtraction, we use the signal-region to sideband ratios obtained 
from Am fits for the data. (Comparison of these ratios for data and simulated BB decays 
shows them to be compatible.) This sideband subtraction correctly accounts for the small 
number of signal decays that populate the Am sideband. 



C. B -> D*iu yields 

Having obtained the distributions in cos9s-D*i of the signal and background components, 
we fit for the yield of D*£u candidates in each w bin. Two representative fits are shown 
for D* + £u in Fig. |TT] and D*°£u in Fig. [T^. The quality of the fits is good, as is agreement 
between the data and fit distributions outside the fitting region. We summarize the observed 
D*lv and D*X£v yields in Fig. p). 

The D*X£u yields correspond to branching fractions of B(B — > D* + X£v) = (0.97 ± 
0.24(stat.))% and B{B -> D*°X£u) = (0.32 ± 0.55(stat.))%. These are somewhat lower 



than past measurements |33|, |34j , but because the analysis is not optimized for these modes 
the systematic uncertainties on these branching fractions are large, of order ±30% for D* + 
and ±60% for D*°, dominated by model dependence in the efficiency to satisfy our lep- 
ton momentum criteria, uncertainty in the correlated and uncorrelated backgrounds, and 
radiative effects in D*£v. 

In order to test the quality of our cos 6b~d*i fits and the modeling of the signal and 
backgrounds, we compare the observed D* energy and lepton momentum spectra with ex- 
pectations for candidates in the signal-rich region \cos 6 b-d*z\ < 1- The cos9b-d*£ fits 



provide the normalizations of the D*£v and D*X£v components. Figure [T3 shows the D 



energy distributions, and Figs. 15 and It show the electron and muon momentum spectra, 
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FIG. 13: The observed (a) and D* + Xl» yields and (b) D* G lv and D*°XlD yields in each 

u; bin. 



repsectively, for D* + £u and D*°£i> candidates. We find good agreement between the data 
and our expectations. 



V. THE \Vcb\ FIT 



The partial width for B — > D*iv decays is given by |35 



as 



dr 

dw 



48tt 3 



(m B — m D *) 2 m 3 Dt Vw 2 — l(w + 1)' 



(6) 



x 



1 + 



Aw 
w + 1 



1 — 2wr + r A 
(1 -rf 



where tub and mo* are the B- and Z)*-meson masses, r = m^./wiB, and the form factor 
•F(u>) is given by 



H 2 n +HI + H1 

4w\ ( l—2wr+r 2 \ 



\i+(^r)( 



h Al (w). 



(?) 



The if j are the helicity form factors and are given by 

w — 1 



1+ , 

1 — r 

\/l — 2wr + r 2 



[1 — R,2(w)] , and 
IT 



w + 1 



Ri(w) 



(8) 
(9) 



The form factor h Al (w) and the form-factor ratios R±(w) = h v (w)/h,A 1 (w) and R-i{w) = 
[hA 3 (vj) + rhA 2 (w)]/hA 1 (w) have been studied both experimentally and theoretically. A 
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CLEO analysis measured these form-factor parameters under the assumptions that 
hA 1 {w) is a linear function of w and that R\ and R 2 are independent of w. CLEO found 



-1 dh 



Ai 



w=l 



= p 

Ri 
R2 



0.91 ±0.15 ±0.06, 



1.18 ±0.30 ±0.12, and 
0.71 ±0.22 ±0.07, 



with the correlation coefficients C(p 2 ,Ri) = 0.60, C(p 2 ,R 2 ) = —0.80, and C(R 1 ,R 2 ) = 
-0.82. 

-Ri(l) and .#2(1) have been computed using QCD sum rules with the results Ri(l) = 1.27 
and i? 2 (l) — 0.8 and estimated errors of 0.1 and 0.2, respectively B7j], in good agreement with 
the (later) experimental results. Ri(w) and R 2 {w) are expected to vary weakly with w. Most 
importantly for this analysis, JF(1)(= /^(l)) is relatively well-known theoretically |15|, 16| , 
thereby allowing us to disentangle it from \ V c b\- 

Recently, dispersion relations have been used to constrain the shapes of the form fac- 
tors [^6|, |38|. Rather than expand the form factor in w, these analyses expand in the 
variable z = (y/w + 1 — y/2) /(y/w + 1 + a/2). The authors of Ref. ^6] obtain 



h Al (w) = h Al (l)[l - 8p 2 z + (53p 2 - 15)z 2 - (231p 2 - 91)^ 3 ], (10) 
R^w) = - 0.12(«; - 1) + 0.05(w - l) 2 , and (11) 

R 2 (w) = R 2 {1) + 0.11(w - 1) - 0.06(w - l) 2 . (12) 

In our analysis, we assume that the form factor has the functional form given in Eqs. |H|- 
12| . We fit our yields as a function of w for |V^|.F(1) and p 2 , keeping R\(l) and R 2 (l) fixed 
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FIG. 15: The electron momentum spectrum for (a) D* + e v candidates and (b) D*°e v candidates 
in the region |cos#b-d*£| < 1 for all w bins combined. 



at their measured values. Our fit minimizes 

10 r a robs V^IO 



\ N obs e ../V-l 2 

x 2 = £- — * . ( 13 ) 



where N° is the yield in the i th w bin, Nj is the number of decays in the j th w bin, and 
the matrix e accounts for the reconstruction efficiency and the smearing in w. 

The efficiency matrix e is calculated using simulated D*iv decays. A matrix element 
represents the fraction of D*lv decays generated in the jth w bin that are reconstructed in 
the ith w bin. To be consistent with our method for finding the cos 9 b-d*£ distribution of 
D*£u decays, described in Section [IV B 61 we subtract the combinatoric background in the 



simulated decays using the Am sideband and the data normalizations. We veto all other 
backgrounds using generator-level knowledge of the simulated events. A single element of 
the efficiency matrix is thus calculated using 

<«, (Sr ■■n,S;""):S j . (14) 

where S^ 9 and Sf lde are the number of non-vetoed candidates reconstructed in the ith w 
bin in the Am signal and sideband regions, respectively, rii is the normalization of the Am 
sideband region, and Sj is the number of D*£u decays generated in the jth w bin. 

The efficiency matrix is nearly diagonal because the resolution in w is about half the bin 
size. The off-diagonal elements are only appreciable for \i — j\ < 1. The resolution becomes 
worse for larger w (see Fig. ^). The efficiency matrix depends not only on the experimental 
selection criteria but also on the form factor. For the cuts described in this paper and using 
the form factor described above, the diagonal elements of the efficiency matrix vary from 
4-14% for D* + £v and from 5-11% for D*°£u. Although we bin in w, the efficiency matrix has 
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FIG. 16: The muon momentum spectrum for (a) D* + \i v candidates and (b) D /i v candidates 
in the region |cos#£-D*£| < 1 for all w bins combined. 



a weak dependence on the slope parameter p 2 . We iterate the fit, reevaluating the efficiency 
matrix for the best-fit value of p 2 . A single iteration is sufficient for convergence. 



In Equation [L3|, Nj is given by 



Nj = 4:f 00 N r ^B D *+B D oT B o 



dw 



dr 

dw 



(15) 



for D* + £i>, where r^o is the B° lifetime [30[|, £>£>.+ is the D* + — > D tt + branching fraction [^0 
Bpo is the D° — > K~n + branching fraction, A^x(45) is the number of T(4S) events in the 



sample, and /oo represents the T(45') — > B B branching fraction. For D 



)*0/ 



dW ~d~w' 



(16) 



where B 
fraction f|30 



D* 



B w o is the 7r c 



77 branching 



is the D*° — > D°7[° branching fraction 

and / + _ represents the T(4S') — > B + B~ branching fraction. The values that 
we use for the B lifetimes and the various branching fractions are listed in Table [VI H For 
B(D° — > A""7r + ), we average the CLEO f39| and ALEPH results after correcting the 
former for final-state radiation (about a 2% correction) to obtain a branching fraction for 
the sum of radiative and non-radiative decays. We exclude the other results included in the 
PDG average because they do not specify their treatment of radiation. 

We first fit D* + £u and D*°£i> separately, allowing as free parameters |V^b|JF(l), p 2 and 

/_i , with the last of these constrained by adding a term (R — Rq) 2 ja 2 R to the x 2 °f Eq. [T3| . 

Here the double ratio R = [/ + _/(l — / + _)](tb+/tbo) is compared to a measurement of the 
same double ratio (_R ± i n Ref. Ell and we have explicitly assumed / o + /+- = 1- The 
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TABLE VII: The lifetimes and the branching fractions used in the \ V c b\ fit. 



t b + (1.653±0.028) ps 

t b o (1.548±0.032) ps 

B(D* + -► D°ir + ) (67.7±0.5)% 
B(D*° -► L> vr°) (61.9±2.9)% 

B(D° -► ^-vr+) + -► (3.89±0.11)% 
£(tt -> 77) (98.798±0.032)% 
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FIG. 17: The results of (a) the D* + £u fit and (b) the D*°£i) fit. The circles and squares are the 
data and the histogram shows the results of the fit, done separately for each mode. 



results of the separate fits are shown in Fig. [17|. For D* + £i>, we find 



= 0.0424 ±0.0018, 
:! - 1.60 ±0.11, and 



P 

x 2 



6.6/8 degrees of freedom (dof.). 



These parameters imply T = 0.0380 ± 0.0019 ps 1 . For D*°£i>, we find 

1^1^(1) = 0.0436 ± 0.0026, 
~ 2 - 1.56 ±0.18, and 



P 

x 2 



9.5/8 dof. 



These parameters imply T = 0.0415 ± 0.0027 ps -1 . The results from D* + £u and D*°£u are 
consistent with each other. 

We also do a combined fit to the D* + £v and D*°£v data. In minimizing, x 2 is the sum of 
the separate D* + £v and D*°£v x 2 's, but including the term constraining (f+-/foo)( T B+/ T B°) 
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only once. The results of the fit are displayed in Fig. 



18 



and the parameter values are: 



|T/ C& |JT(1) = 0.0431 ±0.0013, 
/+_ = 0.521 ±0.012, and 
p 2 = 1.61 ±0.09, with 
X 2 = 16.8/18 dof. 

These parameters give T = 0.0394 ± 0.0012 ps _1 . Not surprisingly, the values of |V^,|.F(1) 
and p 2 are strongly correlated. The correlation coefficients are C(\V c b\J-'(l), p 2 ) = 0.865, 
C(\V cb \T(l),f+-) = 0.130, and C(p 2 ,/+_) = -0.075. 

When we remove the constraint on / H in the fit, we find /+_ = 0.532 ±0.016. This is in 

-0.011 



agreement with the recent CLEO measurement [41], which implies /_| = 0.510±0.017 



-o.oio- 



These results may be compared to our previous analysis |L8| that analyzed a subset 
(approximately 50%) of the current data finding a smaller value for |V^|.F(1). The increase 
may be attributed to several effects. Changes in the measured values of the D° and D* 
branching fractions and B lifetimes cause a 2.3% increase. Inclusion of final-state radiation 
shifts | Vd)|.F(l) by 2.4%. More significantly, use of the improved form factor gives an increase 
of 5.7%. The new form factor has positive curvature, which results in an increase when 
extrapolating to w = 1. The analysis of Caprini et al. []2£j shows there is correlation 
between the curvature and slope of the form factor, making this effect more pronounced for 
the large slope preferred by our data. The remainder of the increase results from the larger 
data sample. 

We test the compatibility of the old and new analyses by restricting the new analysis to the 
same subset of data and fitting using the old form factor. 3 Adjusting for common values for 
D° and D* branching fractions and B lifetimes (Table |V11] ) we find a change in |V^|.F(1) of 
0.0020 ±0.0010 ±0.0022, where the first error is statistical (assessed conservatively assuming 
all candidates in the old analysis are found in the new analysis) and the second error is an 
estimate of the uncorrelated systematic uncertainties. The largest of the latter are due to 
slow pion efficiency, taken to be uncorrelated because of significant differences in the tracking 
algorithms used in the two analyses. We conclude the old and new analyses are compatible 
within the systematic uncertainties. Because our new analysis includes the data reported 
previously and takes advantage of theoretical improvements in the form factor, the results 
reported here supersede our previous results. 

VI. SYSTEMATIC UNCERTAINTIES 

The systematic uncertainties are summarized in Table |V111|. The dominant systematic 



uncertainties arise from our background estimations and from our knowledge of the slow 
pion reconstruction efficiency. 



A. Background uncertainties 

Here we present the systematic uncertainties from our backgrounds. 



Specifically, in |l8) we used a linear form factor h,A x (w) = \ia x (1) [l — p^ii ( w ~~ -0] ; an d Ri = Rz = 1. 
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FIG. 18: The results of the combined fit to the w distribution: (a) the D* + £i> yields (circles) 
with the results of the fit superimposed (histogram) and (b) the D* a li> yields (squares) and fit 
(histogram). In (c) the curve shows the best-fit | V^l^u;), and the circles (squares) are the D* + £v 
(D* ^) yields corrected for efficiency, smearing, and all terms in the differential decay rate apart 
from {VcblTiw). 
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TABLE VIII: The fractional systematic uncertainties, given in percent, for the B° — > D* + £v fit, 
the B — > D*°£u fit, and the combined fit. 





D*+lv fit 




D u 


Iv fit 




Combined fit 




bource 


\Vcb\F(l) 


P 2 


1 




P 


1 




P 2 


1 


Backgrounds 


1.8 


3.0 


1.8 


2.4 


5.0 


2.2 


1.8 


3.1 


1.7 


Reconstruction efficiency 


4.4 


5.0 


A C\ 

4.9 


3.5 


6.2 


6.5 


2.9 


3.2 


4.6 


B momentum & mass 


n 9 


n n 


U.O 


U.O 


U.O 


U.o 


n 1 


n 1 


U.z 


B -> D*XlD model 


0.3 


3.5 


1.2 


1.2 


2.7 


0.5 


0.3 


1.6 


0.9 


Final-state radiation 


0.7 


0.3 


1.1 


0.8 


0.5 


1.2 


0.7 


0.3 


1.1 


Number of BB events 


0.9 


0.0 


1.8 


0.9 


0.0 


1.8 


0.9 


0.0 


1.8 


Subtotal 


4.9 


6.8 


5.8 


4.6 


8.5 


7.3 


3.6 


4.8 


5.4 


tb and branching fractions 


1.5 


0.0 


3.0 


2.8 


0.0 


5.6 


1.8 


0.0 


3.5 


i?i(l) and R 2 (l) 


1.6 


11.7 


1.8 


1.1 


14.3 


1.8 


1.4 


12.0 


1.8 


Subtotal 


2.2 


11.7 


3.5 


3.0 


14.3 


5.9 


2.3 


12.0 


3.9 


Total 


5.3 


13.5 


6.8 


5.5 


16.6 


9.3 


4.3 


13.0 


6.6 



1. Continuum background 

Our estimate of background from e + e _ — > qq is taken from data collected below the 
T(AS). The shape of this background in each w bin is taken from off-resonance data, 
where we scale the energy of the D* and lepton to reflect the difference in the on- and off- 
resonance center-of-mass energies. This scaling applies to computation of w and cos6b-d*£- 
The resulting distribution is scaled by the ratio (C oa / C Q s)(E^ s / E^ n ) , where the first factor is 
the ratio of on- to off- resonance luminosities and the E 2 ratio corrects for the 1/s dependence 
of the hadronic cross section. 

The uncertainty on the normalization is small and has negligible effect on the results 
because the continuum background itself is small. To assess the systematic uncertainty 
from the D* and lepton energy scaling, we compare our results with the scaling to those 
obtained without it. The systematic uncertainties are taken to be half this difference, and 
are 0.03%, 0.2%, and 0.1% for |Kt|.F(l), p 2 , and T(B -> D*£&), respectively. 



2. Combinatoric background 



Our method for combinatoric background subtraction assumes that the cos9b-d*i dis- 
tribution of candidates in the Am sideband matches that of those in the signal region. The 
Monte Carlo simulation should reproduce any differences well since they arise from kine- 
matic effects. We use a sample of 16 million Monte Carlo-simulated inclusive BB events to 
test this assumption. (See also Figs. || and [10] and discussion in Section [IV B 2| .) We per- 
form our analysis on the simulated events twice, once using the combinatoric background 
subtraction procedure outlined in Section IV B 2 and once using the absolutely normalized 
"true" combinatoric background, i.e., D* candidates in the Am signal region that do not 
arise from the decay of a D*, in place of the Am sideband distribution. Use of the "true" 
background instead of the estimate results in a shift in the combined fit of (—1.3 ± 0.9)% in 



1^1^(1), (-2.6 ± 1.3)% in p 2 , and (-0.6 ± 0.9)% in T. Any bias from the use of the Am 
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TABLE IX: The uncertainties due to the combinatoric background. 



Variation \V cb \F{\) {%) p 2 {%) T(B — > D*£v)(%) 



D*+£v fit 

cos 6b-d*£ distribution of sideband candidates 1.6 2.7 1.3 

Am sideband normalizations 0.2 0.3 0.3 

Non-if7r decays (see Table |yD 03 OO 0.6 

Total L6 2?7 L4 

D*°£u fit 

cos 6b-d*£ distribution of sideband candidates 2.2 4.8 1.5 

Am sideband normalizations 0.3 0.6 1.2 

Non-Kn decays(see Table [V]) 0.3 0.0 0.6 

Total 49 2X) 

combined fit 

cos 9b-d*i distribution of sideband candidates 1.6 2.9 1.1 

Am sideband normalizations 0.2 0.3 0.6 

Non-ifvr decays(see Table fVp 0.3 0.0 0.6 

Total L6 2^9 L3 



sideband to estimate the combinatoric background is smaller than the statistical uncertainty 
of the fit to the data. We conservatively assign systematic errors equal to the quadrature 
sum of the shift and its statistical uncertainty, a total of 1.6% for |V^h|jF(l). 

The normalization of the background relies on the fits to the Am distributions. We 
assign an uncertainty for this by repeating our analysis with different functional forms used 
to fit Am. We also include a 0.1% uncertainty because the simulated Am signal peaks are 
shifted a few tenths of an MeV lower than the data. The statistical error on the background 
normalization is included in the statistical errors on our result. 

The final contribution to the systematic uncertainty from our combinatoric background 
estimate comes from the decay modes other than D — > Kit that are reconstructed in our 
m(Kir) signal region. The specific modes were given in Table |IV|. We find the total contribu- 
tion to our D*tv yield from this source is (0.5±0.3)%. We add the yield and its uncertainty 
in quadrature to get a 0.6% uncertainty on our D*tv yield. Because | V C \\T{Y) is proportional 
to the amplitude rather than the rate, its error is half as big. We find the total error on 
|Vd,|.F(l) due to the combinatoric background to be 1.6%. The errors from all components 
are summarized in Table lIX . 



3. Uncorrelated background 

The main source of uncertainty from the uncorrelated background is the normalization of 
the various contributions. Of these, the most important is the normalization of the upper- 
vertex D* decays, which we vary by 50%. Smaller uncertainties arise from the primary and 
secondary lepton rates, the uncertainty in B° — B° mixing, and the uncertainty in the rate 
of exchanging K and tt particles in D*° candidates. The effects of varying these rates are 
summarized in Table [X[ The systematic uncertainties from the uncorrelated background 
estimate are at or below the 1% level. 
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TABLE X: The systematic uncertainties for the uncorrelated background in the separate D* + lv 
and D* a iu fits and the combined fit. The upper-vertex contribution, the lepton normalizations, 
and Xd are treated as completely correlated between the two modes; all others are uncorrelated. 



Variation 


\V cb \f(l)(%) 


p 2 (%) 


T(B -» D*£u)(%) 


D*+lu fit 








Upper-vertex D* 


0.7 


0.9 


0.4 


Other 


0.2 


0.2 


0.3 


Total 


0.7 


1.0 


0.5 


D*°iv fit 








Upper-vertex D* 


0.6 


0.9 


0.4 


Other 


0.5 


0.8 


0.3 


Total 


0.8 


1.2 


0.5 


Combined fit 








Upper-vertex D* 


0.6 


0.9 


0.4 


Other 


0.3 


0.3 


0.3 


Total 


0.7 


1.0 


0.5 



4- Correlated background 

We assess the uncertainty arising from the correlated background by varying the branch- 
ing fractions of the contributing modes simultaneously by 50%. Since this is a small back- 
ground, this variation has little effect on |Vd,|.F(l), and the uncertainties are 0.1%, 0.6%, 
and 0.8% on p 2 , and T(B -> D*£v), respectively. 

5. Fake lepton background 

We vary the measured electron and muon fake rates separately by 50%. This is conser- 
vative, but it has also almost no effect on our result; the total uncertainty on |Vd>|.F(l) is 
0.02%, while the uncertainties on p 2 and T(B — > D*£u) are 0.3% and 0.2%, respectively. 

B. Slow 7r reconstruction uncertainty 

The largest source of uncertainty for the analysis is the efficiency for reconstructing the 
slow pion from the D* decay. Because of the small energy release in D* decays, the daughter 
pion has low momentum and travels approximately in the direction of the parent D*. For 
our signal decays, the momentum range of the slow pion is to about 250 MeV/c. Note 
also that w = Ed* /mo* hi the B rest frame, so the slow-pion momentum is correlated with 
w [see Fig. 0(a)]. 

Charged and neutral slow-pion reconstruction efficiencies depend very differently on w. 
Charged pions with momenta less than 50 MeV/c do not penetrate far enough into the 
tracking chamber to be reconstructed; the slow-pion reconstruction efficiency is therefore 
low near w — 1 and increases rapidly over the next few w bins as the pion momentum 
increases. Neutral slow pions, on the other hand, decay to two low-energy photons (30- 
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FIG. 19: (a) The correlation between the slow pion momentum and w; (b) The reconstruction 
efficiency as a function of w for charged slow pions (solid circles) and neutral slow pions (open 
squares) . 

230 MeV). The lowest-momentum 7r°'s decay almost back-to-back, depositing about equal 
energy in the calorimeter. As the 7r° momentum increases, the Lorentz boost pushes some 
of the photons below our minimum energy requirement of 30 MeV. The neutral slow-pion 
efficiency therefore drops slowly as w increases. The slow-pion efficiencies for both charged 
and neutral 7r's from B — > D*iv decays are shown as a function of w in Fig. |l9l(b). 

Because we rely on Monte Carlo simulation to estimate the slow-pion efficiencies, we in- 
vestigate possible differences between the simulation and performance of the CLEO detector 
in order to estimate the systematic uncertainty of slow-pion reconstruction. We consider 
the effect of nearby tracks and showers on slow-pion reconstruction, comparing the efficiency 
of data and Monte Carlo-simulated BB events to limit a systematic error due to a differ- 
ence in the "event environment" in data and simulated events. We also consider how much 
imperfect knowledge of detector material can affect reconstruction efficiency through pion 
range-out, multiple scattering, hadronic interaction, or photon conversions. Finally, we vary 
parameters of the detector simulation for the drift chambers and calorimeter to estimate the 
contribution to the systematic uncertainty. 

Monte Carlo-simulated events may have a different number of drift chamber hits or 
calorimeter showers than the data. The detector activity (track fragments or showers) 
near a candidate slow pion can affect the reconstruction efficiency. To evaluate the impact 
of these "environment effects" on the slow-pion reconstruction efficiency, we insert Monte 
Carlo-generated slow-pion tracks or showers with kinematic distributions appropriate to 
D*£u decay into samples of hadronic events selected from our data and from simulated BB 
events. In each w bin, we compare the reconstruction efficiency for the tracks embedded into 
data and simulated events. For D* + £u, the efficiency difference is small, Ae/e = (0.2±1.1)% 
integrated over w. Likewise, the effect of event environment is small for D*°£u, where we 
find a net efficiency difference Ae/e = (—0.6 ± 1.1)%. The uncertainties here are from the 



29 



statistics of the data and Monte Carlo comparison. We measure the impact of the event 
environment by using the measured data-Monte Carlo efficiency difference in each w bin to 
modify the efficiency matrix in Eq. [13] and repeating the fit. The slow-pion efficiency may 
depend on the track or shower multiplicity, which is increased by one or two, respectively, 
by the embedding study; we find no statistically significant evidence of this in our studies, 
but we include a small uncertainty [0.3% on (^(^(l)] to cover this effect. 

To estimate the uncertainty due to our imperfect knowledge of the detector material 
inside the outer boundary of the tracking chambers, we vary the material description of 
the detector by 10% in our simulation and remeasure the slow-pion efficiencies. This 10% 
variation of material is based on a study that compared the polar angle distribution of 
e + e~ — > 77 events in data and simulation. We then repeat the \V c b\ fit using these new 
efficiencies and take the excursions of |V^b|jF(l), p 2 , and T(B — > D*lv) as the uncertainty. 

In a similar way, we estimate the uncertainty due to our tracking chamber and crystal 
calorimeter simulation. For charged slow pions, performance of the tracking devices is 
essential. Differences in hit efficiency and single-hit resolution between data and Monte 
Carlo simulation can result in a difference in measured efficiency. The tracking simulation 
parameters are tuned using an independent sample of charged tracks. We vary the tracking 
chamber hit resolutions by amounts determined from residual distributions in these data, 
and we vary hit efficiencies according to observed differences in the data and simulated hit 
efficiencies. 

For neutral slow pions, performance of the calorimeter is important. Here we consider 
differences in the ^(77) and transverse shower profile distributions used for 7r° reconstruc- 
tion. We calibrate the calorimeter energy scale at high energy (1-5 GeV) using showers from 
QED event samples (e + e~ — > e + e~ , e + e~ — > 77, and e + e~ — > e + e~ r y). We check this scale 
with a sample of 7r° — > 77 and 77 — > 77 candidates, which should peak at the known 7r° and 
77 masses. For low-energy showers there can be residual gain mismatches from nonlinearities 
and noise. Accordingly, we adjust the calorimeter noise and dispersion of crystal gains in the 
simulation so that it reproduces the transverse spatial distributions of the photon showers 
and the ^(77) distribution for an independent sample of low-momentum 7r°'s. We vary 
the noise and gain dispersion parameters in the simulation within a range determined from 
the data to assess the systematic uncertainty. Photons that convert and begin to shower 
just in front of the calorimeter will have degraded resolution. We vary the material de- 
scription between the outer tracking chamber boundary and the calorimeter crystals by a 
conservative 15% to determine its contribution to the uncertainty for slow-7r° reconstruc- 
tion. The transverse spatial extent of photon showers varies with the low-energy cutoff in 
the shower simulation. To assess the uncertainty from the cutoff in our simulation, we lower 
the minimum energy for photon simulation by a factor of 10, from 1 MeV to 100 keV. 

Finally, we assess the systematic uncertainty due to requiring a D* + and D° vertex in 
the D* + lv analysis by performing the analysis without vertexing. In the separate D* + £u 
fit without vertexing, the result for |V^h|jF(l) shifts by (2.0 ± 1.8)%, where the uncertainty 
takes into account correlations between analyses with and without vertexing. We take the 
quadrature sum of the shift and its uncertainty as a systematic error. 

We find that the largest contributions to the uncertainty on |V^f,|jF(l) come from the 
material description (1.3%), the effects of vertexing (1.5%), and the minimum energy for 
photon simulation (0.6%). The statistical uncertainty from data and Monte Carlo compar- 
isons also contributes (1.3%). The given uncertainties apply to the combined fit. Table [XI 
summarizes the uncertainties on slow-pion reconstruction. 
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TABLE XI: The systematic errors from the slow-pion reconstruction efficiency for the separate 
D* + lv and D*®li> fits and the combined fit. We take the uncertainty from the material description 
to be correlated between the two modes; all other errors are uncorrelated. 

Mode |V c6 |^(l)(%) p 2 {%) T(B -> D*lv){%) 



D*+£v fit 

Material description 2.6 3.2 2.1 

Tracking chamber hit efficiency 0.6 0.2 1.4 

Vertexing 2.7 2.7 2.9 

Other uncertainties 0.8 1.1 0.7 

Statistics (environment) 1.7 2.5 1.4 

Total 42 Ko IT 
D*°lv fit 

Material description 1.1 3.0 0.6 

Photon cutoff 1.5 0.9 2.3 

Other uncertainties 1.2 2.9 5.0 

Statistics (environment) 2.1 3.3 2.7 

Total 3T Kl 6^2 
Combined fit 

Material description 1.3 1.5 1.2 

Tracking chamber hit efficiency {D* + £,u only) 0.3 0.2 0.9 

Vertexing (D*+£u only) 1.5 1.6 1.7 

Photon cutoff (D*°£i> only) 0.6 0.2 0.9 

Other uncertainties 0.9 1.0 1.8 

Statistics (environment) 1.3 2.1 1.3 

Total 2J3 3A 3~T 



C. Sensitivity to i?i(l) and R 2 (l) 

The form factor ratios -Ri(l) and -R 2 (l) affect the lepton spectrum and therefore the 
fraction of decays satisfying our 0.8 GeV/c electron and 1.4 GeV/c muon momentum re- 
quirements. They also affect the relative contributions of the three D*£u form factors, and 
therefore can affect the form-factor slope p 2 . 

To estimate the uncertainty due to the measurement errors on Ri(l) and -R 2 (l), we use 

2 A OP dP 

ffp = ^il % (17) 

where P stands for the parameter [|Vd,|.F(l), p 2 , or T(B — > D*iv)\ whose uncertainty we 
are calculating, En = af and = pij<Ti<Tj, where p i2 = —0.82 is the correlation coefficient 
from the Ri(l) and -R 2 (l) measurement |^6[. We compute the partial derivatives dP/dRiil) 
by shifting Ri and repeating our analysis. We find an uncertainty on |V c fe|jF(l) from this 
source of 1.4%, and a substantial uncertainty on p 2 of 12%. 
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TABLE XII: The results from separate analyses using only D*e v or D* ji v. The errors are 



statistical only. In these fits /_| is constrained using Ref. ^l|. 



Mode 


it/ i -7-/1 "\ t<y \ 


P (/o) 


1 {rS — > JJ lu){/o) 


JJ ^iv 








D* + e~v 


0.0420±0.0023 


1.65±0.14 


0.0363±0.0021 




0.0448±0.0026 


1.69±0.15 


0.0404±0.0025 


D*°ii> 










0.0409±0.0032 


1.41±0.24 


0.0396±0.0030 




0.0474±0.0040 


1.80±0.26 


0.0423±0.0042 


combined 








D*e~9 


0.0420±0.0018 


1.60±0.12 


0.0374±0.0015 


D*pTv 


0.0457±0.0021 


1.73±0.13 


0.0411±0.0019 



D. Other uncertainties 

We considered the following minor sources of systematic uncertainty, summarized in 
Table [Vm]. 

The efficiency for identifying electrons has been evaluated using radiative Bhabha events 
embedded in hadronic events, and has an uncertainty of 2.6%. Similarly, the muon identifi- 
cation efficiency has been evaluated using radiative mu-pair events, and has an uncertainty of 
1.6%. We determine the total uncertainty from lepton identification by adding in quadrature 
the shift in results when repeating the analysis with electron and muon efficiencies varied 
by their momentum-dependent uncertainties. Separate electron and muon analyses of our 



data give the results shown in Table |XII| . Including the systematic uncertainties on lepton 
identification, the separate electron and muon results are consistent at the 35% confidence 
level. 

The B momentum is measured directly in the data using fully reconstructed hadronic 
decays, and is known on average with a precision of 0.0016 GeV/c. Variation of the mo- 
mentum in our reconstruction slightly alters the cos6b-d*£ distribution that we expect for 
our signal, and it therefore changes the yields obtained from the cos 6 b-d*£ fits. Likewise, 
CLEO has measured the B°- and S + -meson masses [|4] and when we vary them within 
their measurement errors, we find a small effect on the yields. 

We determine the tracking efficiency uncertainties for the lepton and the K and 7r forming 
the D° in the same study used for the slow pion from the D* + decay. These uncertainties 
are confirmed in a study of 1-prong versus 3-prong r decay events from our data sample. 

The final-state radiation model has a small effect on our D*iv yields because it affects 
the D*lv cos 9 B _ D *i distributions. Because we require p e > 0.8 GeV/c and > 1.4 GeV/c, 
the model also affects the D*£u efficiency. The final-state radiation model is estimated by 
the authors of Photos to be accurate within 30% P7|. We determine our sensitivity to 



the model by repeating our analysis without including radiative D*lv decays in our D*tv 
Monte Carlo. We then take 30% of the change to our results as our uncertainty. 

Finally, our analysis requires that we know the cosQb-d*i distribution of the D*X£u con- 
tribution. This distribution in turn depends on both the branching fractions of contributing 
modes and on their form factors. Variation of all of these branching fractions and form 
factors is not only cumbersome, but also out of reach given the poor current knowledge of 
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FIG. 20: The cos6b-d*£ distribution of B —> D* n£ij (solid histogram) and B — * D\lv (dashed 
histogram) decays contributing to the D*Xlv sample for D*°lu. The histograms are normalized 
to equal area. 

these modes. Instead, we note that the B — > D*n£u and B — > D\£v modes are the ones with 
the most extreme cos8b-d*i distributions (the largest mean and the smallest). These distri- 
butions are shown in Fig. |20|. We therefore repeat the analysis, first using only B — > D*itiv 
to describe our D*Xtv decays and then using only B — > D\tv to describe these decays; we 
take the larger of the two excursions as our systematic error. 



VII. CONCLUSIONS 

We have fit the w distribution of B — > D*iv decays for the slope of the form factor and 
1^1^(1). For the combined D* + £u and D*°£v fit, we find 



Including the systematic uncertainties we compute a correlation coefficient 
C(\V c b\J I '(l), p 2 ) = 0.22. Figure ^T] shows the total error ellipse for this measurement. The 
best-fit parameters imply the decay rate 



We recover the branching fractions from the rate by dividing by the appropriate I?-meson 
lifetimes. These results are sensitive only to the ratio of B + to B° lifetimes. They are 



V cb \F{l) = 0.0431 ±0.0013 ±0.0018, and 
p 2 = 1.61 ±0.09 ±0.21. 



T = 0.0394 ± 0.0012 ± 0.0026 ps 



-l 



B(B°^D* + £u) = (6.09 ±0.19 ±0.40)% and 
B(B~ ->• D*°£u) = (6.50 ±0.20 ±0.44)%, 
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FIG. 21: The error ellipse for the combined D* + lv and D* G £v measurement, including statistical 
and systematic uncertainties. 



where the errors are completely correlated. 

A recent lattice calculation yields |42|] 
of +0.007. This value is consistent with .F(l) 



= 0.9191 

ao35 a ft er a Pply m g a QED correction 
0.913 ± 0.042, the evaluation of the authors 



of Ref. [43], but is more precise. Using the lattice value of .F(l), our result implies 

0.0469 ± 0.0014(stat.) ± 0.0020(syst.) ± 0.0018(theo.). 



V, 



cb\ 



Since full radiative corrections have yet to be calculated for ^-"(1), it is ambiguous how best 
to treat radiative decays in the analysis. We include radiative decays in our signal. 
This value of 



decays 



15 



\Vcb\ is consistent with previous values obtained from D 
but is somewhat higher. However, we note our ability to reconstruct cos Qb-d*i 
makes our analysis approximately four times less sensitive to the poorly known D*X£u 
background, and furthermore allows us to constrain it with the data. This value of \V c b\ 
is also somewhat higher than that obtained using inclusive semileptonic B decays |4"8|| . If 
confirmed, this discrepancy could signal a violation of quark-hadron duality. A larger value 
of \V c i\ shifts constraints on the CKM unitarity triangle from \ V u b/V c b\ and CP violation in 
the neutral kaon system, somewhat reducing expectations for indirect CP violation in the 
B system. 
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